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Measurements of J/V^ production by the PHENIX experiment in p+p, d+Au, Cu+Cu and 
Au+Au collisions at y^Siviv=200 GeV are reviewed. The results show a suppression beyond 
what can be explained by cold nuclear matter effects in the most central Au+Au and to a 
lesser extent in Cu+Cu collisions. In addition, the suppression observed at mid rapidity in 
Au+Au is smaller than at forward rapidity, a tendency opposite to what is expected from the 
higher energy density at mid rapidity. Regeneration, a possible explanation, can be tested by 
measuring the elliptic flow parameter V2 of J/ip. 



1 Introduction 



J/t/j suppression is considered to be one of the key probes of the Quark Gluon Plasma (QGP) 
formation in heavy ion collisions. Color screening was proposed'!] as a mechanism leading to 
anomalous suppression beyond normal hadronic absorption if J/'0s are created in a deconfined 
medium. The CERN SPS experiments NA38, NA50 and NA60 were the first to investigate this 
phenomenon by measuring J ftp suppression in a variety of colliding systems and energies. The 
results show a statistically significant anomalous suppression in central Pb+Pb ^ and In+In^ 
collisions, that can be interpreted in terms of melting in the QGP. 

The PHENIX experiment at RHIC has also measured the production of J/'0 in a variety of 
colliding systems, and provided further insights by exploring this phenomenon at higher energies. 

are detected in PHENIX through their dielectron decay at mid rapidity {\y\ < 0.35) 
and through their dimuon decay at forward rapidity (1.2 < \y\ < 2.4). J/'0 suppression is 
characterized by a ratio called the nuclear modification factor, obtained by normalizing the J/i/j 
yields in heavy ion collisions {dNAs) by the J/'^/^ yields in p+p collisions at the same energy 




(dNpp) times the average number of binary inelastic nucleon-nucleon collisions (< Ncoii >)'• 

o / N dNAB{y^PT)/dydpT ... 

< Ncoii > dNpp{y,pT)/dydpT 

If the heavy ion collision is a superposition of independent ^coll inelastic nucleon-nucleon colli- 
sions, Kab will be equal to unity, whereas it will be larger than one in case of enhancement and 
lower than one in case of suppression. 

At typical RHIC energies {^/snn — 19.6 - 200 GeV), J/t(js are dominantly produced through 
gluon fusion. The J/^p yield is therefore sensitive to gluon shadowing §. Part of the ground state 
charmonia yield also comes from feed down of excited states (^^ and Xc)? ^.nd can contribute up 
to ^ 40% to the total J/^ yield. In subsequent stages of the collision involving heavy ions, there 
are a number of competing mechanisms that can enhance or suppress the J/tp yield. The two 
major contributors to the suppression are absorption by nuclear fragments from incident nuclei, 
and an eventual melting in the QGP. Finally it is not impossible that a pair of uncorrelated c 
and c quarks that are close enough in phase space recombine to form a bound charmonium state 
and enhance the J/tp yield. 



2 Baseline and cold nuclear matter effect measurements 



The differential cross section of J/ip in p+p collisions as a function of rapidity measured by 
PHENIX is shown in Fig. [T] (left) In addition to providing normalization cross sections 
essential for the calculation of Kaa as in Eq. [H J/'0 measurements in p+p collisions constrain 
the poorly understood J ftp production mechanism. In the same figure, predictions from various 
LO and NLO calculations are shown. The current precision does not discriminate between the 
models, but there is potential for improvement through better precision in cross section and 
additional information from J ftp polarization measurements. 

Nuclear absorption and shadowing, collectively referred to as cold nuclear matter effects 
( CNM) can be constrained by measurements in proton (or light ion) on heavy ion collisions. In 
PHENIX this was performed in deuteron-gold collisions. The resulting suppression ratio ^dAu 
is shown in Fig. [T] (right) ^ as a function of rapidity, where the positive rapidity coincides with 
the deuteron going direction. J/^s detected in different rapidity regions probe specific gluon 
X2 regions |. Forward rapidity corresponds to X2 ^ 0.002 - 0.01 where the depletion due to 
shadowing is important whereas backward rapidity corresponds to X2 ^ 0.05 - 0.2 where a slight 
enhancement due to anti-shadowing is expected. The rapidity dependence of R^^^ therefore 



^ Shadowing refers to the depletion of low momentum partons in nucleons bound in nuclei as compared to free 
nucleons. 

^ By X2, we refer to the parton longitudinal momentum fraction in the nucleus. 



reflects the gluon shadowing, whereas the global vertical scale is determined by the amount of 
normal absorption. 

To quantitatively disentangle the shadowing component from the absorption component, a 
rapidity dependence using two shadowing schemes, EKS and NDSG was fitted to R^^^ leaving 
the overall vertical scale a free parameter to account for the absorption^. J ftp absorption cross 
sections of 2.8til ^nd 2.2^^;^ mb were obtained for EKS and NDSG schemes respectively. 
This is in agreement with the absorption cross section reported by the SPS of 4.2 ± 0.5 mb^ 
but such a comparison should not be taken at face value because shadowing is not taken into 
account in the SPS | absorption cross sections evaluation. 



3 Anomalous suppression in heavy ion systems 

PHENIX has also measured J/tp suppression in Au+Au^ and Cu+Cu"'-^ colhsions at ^/sNN = 
200 GeV. The J/^p Raa in Au+Au colhsions as a function of the number of participants N^^^^, 
at forward and mid rapidity ranges is shown in Fig. [2] (left) ^ together with data points from 
NA38, NA50 and NA60 experiments. The R^^ goes down to ^ 0.2 for the most central Au+Au 
collisions (large ^part)i approaches unity for peripheral ones (small ^part)- To see the extent 
of anomalous suppression, extrapolations of the CNM and shadowing constraints obtained from 
d+Au measurements were calculated using a model dependent method which assumes the above 
mentioned shadowing schemes as well as with a data driven method which has minimal model 
dependence. The result-^ from both methods is a statistically significant suppression beyond 
CNM extrapolation in the most central forward rapidity Au+Au collisions, less pronounced at 
mid rapidity Au+Au or in Cu+Cu collisions. 
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Figure 2: 3/ip Raa vs. Npart at SPS compared to RHIC (left). J/ip V2 vs. pr measurement by PHENIX (right). 

The data show two features that contradict local density induced suppression models. The 
mid rapidity suppression is lower than the forward rapidity suppression (cf. Fig. [2] (left)), despite 
experimental evidence § that energy density is higher at mid rapidity than at forward rapidity. 
The same remark holds for the comparison between Kaa at mid rapidity in PHENIX and Raa 
at SPS I (cf. Fig. [2] (left)). The two are in agreement within error bars, a surprising result 
considering that the energy density reached at RHIC is larger than the one reached at SPS. A 
number of explanations have been put forth, including sequential melting, where only '0' an d Yr 
are dissociated leading to a suppression of only the feed down component of the J/^ yieldlJ^, 
and gluon saturation that leads to a lower charm quark yield at forward rapidity^. 



^ Taking into account nuclear PDF modification would increase the SPS absorption cross section, because the 
SPS rapidity corresponds to the anti-shadowing regime, requiring more absorption to account for the observed 
suppression. 

^ The rapidity density of charged particles which increases with the deposited energy peaks at mid rapidity^^. 
^ Care must be taken when comparing with SPS, because the CNM effects are not the same at the two energies. 



4 Regeneration 



A strong regeneration of J/'0 from uncorrelated c and c quarks is another good candidate to 
explain the tendency of Raa as a function of rapidity at RHIC. This is supported by the high 
charm quark yield measurements (^10 cc pairs are created in the most central Au+Au colli- 
sions). A number of model predictions that incorporate regeneration have been proposed^ and 
all of them reproduce qualitatively the rapidity dependence of J/'0 Raa observed by PHENIX. 

However, important inputs to regeneration models such as the precise number of cc pairs 
available for recombination and the phase space conditions for recombination to take place are 
poorly constrained. It is thus very compelling to have a direct experimental check of regenera- 
tion. The J/'0 elliptic flow is one candidate. Elliptic flow refers to the azimuthal angle correlation 
of particle emission with respect to the reaction plane orientation 0. It is quantified by the second 
Fourier coefficient V2 of the azimuthal angle distribution of identified particles. The measured 
V2 of electrons from D and B meson decays is remarkably high ^. This is believed to originate 
from the elliptic flow of underlying charm and beauty quarks. J/'^s from recombination should 
inherit the charm quark flow, resulting in a higher V2 than the case of direct production in hard 
collisions. 

The first measurement of J/'0 V2 at RHIC energy was performed by PHENIX at mid rapidity 
and is shown in Fig. [2] (right) as a function of transverse momentum. Predictions from models 
that assume various amounts of recombination from none to full coalescence at freeze out are 
plotted together. Data points are compatible within the error bars simultaneously with zero flow 
as well as with the model that predicts maximum flow. This result should therefore be seen as 
proof of principle of the feasibility of J/'0 V2 measurements. There is still room for improvement 
using already existing data, but it is to be noted that a much larger sample will probably be 
needed to be able to distinguish between the different models. 
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The reaction plane is the plane defined by the beam axis and the line joining the center of colliding nuclei. 
It is measured in PHENIX from azimuthal angle distribution of charged particles close to beam rapidity. 



